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ABSTRACT 

Numerical MHD simulations play increasingly important role for understand¬ 
ing mechanisms of stellar magnetism. We present simulations of convection and 
dynamos in density-stratified rotating spherical fluid shells. We employ a new 
3D simulation code for the solution of a physically consistent anelastic model of 
the process with a minimum number of parameters. The reported dynamo simu¬ 
lations extend into a “buoyancy-dominated” regime where the buoyancy forcing 
is dominant while the Coriolis force is no longer balanced by pressure gradients 
and strong anti-solar differential rotation develops as a result. We find that the 
self-generated magnetic fields, despite being relatively weak, are able to reverse 
the direction of differential rotation from anti-solar to solar-like. We also hnd 
that convection flows in this regime are significantly stronger in the polar re¬ 
gions than in the equatorial region, leading to non-oscillatory dipole-dominated 
dynamo solutions, and to concentration of magnetic field in the polar regions. 
We observe that convection has different morphology in the inner and at the 
outer part of the convection zone simultaneously such that organized geostrophic 
convection columns are hidden below a near-surface layer of well-mixed highly- 
chaotic convection. While we focus the attention on the buoyancy-dominated 
regime, we also demonstrate that conical differential rotation profiles and persis¬ 
tent regular dynamo oscillations can be obtained in the parameter space of the 
rotation-dominated regime even within this minimal model. 

Subject headings: dynamo, convection, magnetohydrodynamics (MHD), stars: 
magnetic field. Sun: magnetic fields. Sun: rotation 
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1. Introduction 

Recent progress in observational and compntational capabilities have led to snbstantial 
advances in onr nnderstanding of the origins of stellar and planetary magnetism and 
variability. The key role is played by dynamo processes driven by convection in the presence 
of rotation. This process is particnlarly complex in the case of tnrbulent convection in 
stellar and planetary envelopes. The interaction of convection, magnetic held and rotation 
resnlts in variations of the rotation rate with depth and latitnde, called differential rotation, 
and also in large-scale meridional hows. The diherential rotation is a crucial part of dynamo 
mechanisms. It has been measured quite accurately for the solar surface by tracking motion 
of various features and also through analysis of the Doppler shift of spectral lines. Moreover, 
helioseismology data from the SOHO and SDO space missions and from the ground-based 
network GONG have provided measurements of the internal rotation (e.g. Schou et al. 
1998) and of the meridional circulation of the Sun (Zhao et al. 2013). 

Recently, high-precision spectroscopic and photometric observations have enabled 
measurements of the diherential rotation on other stars. The most prominent feature of 
the solar diherential rotation is that the equatorial zone rotates faster than higher latitude 
regions. For other stars, this type of diherential rotation is called solar-like rotation. When 
the equator rotates slower than the polar regions then such rotational prohle is called 
anti-solar. Anti-solar diherential rotation has been observed by Doppler imaging techniques 
on several K-giant stars (Strassmeier et al. 2003; Kovari et al. 2015). Recent analysis of 
the high-precision light curves from the Kepler mission for a sample of 50 G-type stars by 
Reinhold & Arlt (2015) found 21-34 stars with the solar-like diherential rotation and 5-10 
stars with the anti-solar rotation. The latter work awaits conhrmation from independent 
studies. 

One of the hrst theories of non-linear convection in rotating shells developed in the 
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Boussinesq approximation by Busse (1970), also (Busse 1973), demonstrated that the 
dynamical effects of rotation on convection are a primary mechanism of stellar differential 
rotation. This theory revealed that the basic properties of differential rotation primarily 
depend on the supercritical value of the Rayleigh number which measures the strength of 
the buoyancy forces and determines the magnitude of convection motions. For a relatively 
small supercritical Rayleigh numbers convection develops mostly in the equatorial region 
in the form of convective rolls (also known as ’’banana cells”) oriented along the rotation 
axis. Angular momentum transport by these cells causes solar-like differential rotation. 
This regime of weakly supercritical Rayleigh numbers corresponds to small Rossby numbers 
dehned as the ratio of the rms convective velocity to the mean rotational velocity. This 
regime is called rotationally-dominated. Early numerical simulations in the Boussinesq 
approximation by Gilman (1976) conhrmed these results, and also found that at high 
supercritical Rayleigh numbers when convection develops at all latitudes differential rotation 
may become anti-solar. This regime is characterized by large (typically greater than 1) 
Rossby numbers and is called buoyancy-dominated.. While the rotationally-dominated 
regime is likely to occur in the deep convection zone where convective velocities are 
small, Gilman & Foukal (1979) noticed that on the Sun a typical supergranulation 
turn-over time is much shorter than the rotation period, and thus the convection is in 
the bouyancy-dominated regime. This suggests that on the Sun a combination of the two 
regimes takes place, resulting in a complicated differential rotation prohle. 

Indeed, the differential rotation prohle of the Sun determined by helioseismology (e.g. 
Schou et ah 1998) turned out to be quite different from the theoretical predictions. It is 
characterised by nearly radial orientation of angular velocity at midlatitudes (the so-called 
“conical prohle”). The angular velocity also decreases monotonically from the equator 
to the poles by about 30%. In addition, helioseismology inferences reveal two narrow 
rotational shear layers at the boundaries of the convection zone: the so-called tachocline at 
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the bottom and a near-surface shear layer at the top. 

Substantial efforts have been made to reproduce the observed solar rotation in 
three-dimensional numerical simulations. These efforts are reviewed by Miesch (2005) 
and more recently by Brun et ah (2014). Historically, simulations have been successful in 
reproducing the decrease of angular velocity from equator to poles. However, it is a major 
discrepancy that most simulations tended to feature angular velocity contours parallel to 
the rotation axis (the so-called “cylindrical prohles”), rather than conical prohles. The 
prohle of the differential rotation is governed by the zonal component of the vorticity 
equation and its analysis suggests that cylindrical prohles may be avoided if (a) baroclinic 
forcing with non-vanishing latitudinal entropy gradient exists, or if (b) sufhciently strong 
Reynolds stresses or (c) Lorentz forces develop (e.g. Miesch 2005; Miesch et ah 2006). 
Pursuing alternative (a), Miesch et al. (2006) imposed a latitudinal gradient of entropy as 
a bottom boundary condition in global convection simulations. They were so successful in 
hnding conical differential rotation prohles that many groups have subsequently adopted 
this model as standard (e.g. Browning et al. 2006; Nelson et al. 2013; Fan & Fang 2014). 
While there is physical justihcation for a non-vanishing latitudinal entropy gradient, this 
model assumption must be regarded with caution as the imposed latitudinal variations 
in entropy have not been generated self-consistently by the convective motions in the 
simulations. Large amplitude Reynolds stresses required to realize alternative (b) can be 
achieved by driving convection more strongly. Because most early studies argued that the 
solar magnetic held cannot exert a substantial control over diherential rotation, alternative 
(b) has been explored most often. A recent work in this direction is the attempt of Guerrero 
et al. (2013b) to hnd a regime close to the real solar rotation by varying the gradient of 
the background specihc entropy and the frame rotation rate in a numerical anelastic model. 
These variations correspond to an increasing ehective Rayleigh number and thus to a more 
strongly driven convection. However, the simulations did not hnd an intermediate regime 
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with a “conical” profile. Instead they surprisingly showed that the transition between the 
solar and anti-solar rotational prohles is rather sharp. Most recently similar results have 
been reported by Gastine et al. (2013, 2014); Kapyla et ah (2014); Mabuchi et ah (2015), 
and Karak et al. (2015). 

Alternative (c), namely that dynamo effects are crucial in shaping the prohle of 
differential rotation, is much less explored. A signihcant number of dynamo simulations 
have been published previously, but most of them aimed to reproduce and explain solar 
magnetic features and activity cycles, (e.g. Browning et al. 2006; Ghizaru et al. 2010; Nelson 
et al. 2013). Only recently few studies, of which we mention (Fan & Fang 2014; Mabuchi 
et al. 2015), have appeared that aim to investigate whether the solar magnetic held plays an 
active role in shaping the solar differential rotation prohle. These studies report evidence in 
support of this hypothesis. This is hardly surprising as it is well-established that the main 
ehect of self-sustained magnetic held on convection is to suppress diherential rotation (e.g. 
Grote & Busse 2001; Busse et al. 2003; Simitev & Busse 2005). This hnding is conhrmed 
by Aubert (2005) and Yadav et al. (2013) who also provide scaling laws for this ehect. 
However, the ehects of dynamo action on convection in the buoyancy-dominated regime are 
less well explored. Studies include the afore mentioned papers by (Fan & Fang 2014; Karak 
et al. 2015; Mabuchi et al. 2015) who report a number of similar results but their models 
are rather diherent from each other and in all cases also include factors that contribute to 
baroclinic forcing such that hypothesis (c) is not addressed in isolation. 

In this context our paper has several goals. We wish to study the ehects of self-generated 
magnetic helds on the convective hows of a density-stratihed huid in rotating spherical 
shell using a minimal self-consistent model. We wish to focus the attention on dynamos 
near the transition to buoyancy-dominated convection but we also report selected results 
in the rotation-dominated regime. To this end, we employ the so-called Lantz-Braginsky 


7 


anelastic approximation (Braginsky & Roberts 1995; Lantz & Fan 1999; Jones et al. 2011). 
It has the advantages that the dynamics of the system depends on a minimal number 
of non-dimensional parameters while ad-hoc parametrizations of physical effects that 
may be used to better “£t” observations are excluded. We present a newly implemented 
numerical simulation code for the solution of the Lantz-Braginsky anelastic equations, as 
well as code validation results based on published benchmark solutions against four other 
independently-developed codes (Jones et al. 2011). Since the dynamo processes operating 
in the Sun and stars are still subjects to controversies with various competing proposals, we 
believe that it is important that our results are readily reproducible by other groups. 

In Section 2 we introduce the mathematical model based on the anelastic approximation, 
and discuss the numerical method and diagnostic output. In Section 2.5 we present the 
benchmark validation results. In Section 3 we discuss the properties of dynamos in the 
buoyancy-dominated regime of convection. In Section 4 we demonstrate that conical 
differential rotation prohles and persistent regular dynamo oscillations can be obtained in 
the rotation-dominated regime. In Section 5, we present a summary of our main results and 
compare them to related recent studies. We also outline topics for future research. 


2. Mathematical model and numerical method 

We consider an electrically conducting, perfect gas conhned to a spherical shell. The 
shell rotates with a hxed angular velocity Qk about the vertical axis and an entropy 
contrast AS* is imposed between its inner and outer surfaces. 


2.1. Anelastic governing equations 


Assuming a gravity field proportional to l/r^, a hydrostatic polytropic reference state 
exists of the form 


p = T = T,C, P = PcC^\ C = Co + Cid/r, (1) 

with parameters Cq = (2Co - rj - 1)/{1 - rj), ci = (1 + r]){l - Co)/(l - vf, Co = 

{jj + l)/(? 7 exp(Ap/n) + 1). The parameters Pc, Pc and Tc are reference values of density, 
pressure and temperature at the middle of the shell. The gas polytropic index n, the 
density scale height number, Np, and the radius ratio, p, are dehned below. Convection and 
magnetic held generation are described by the equations of continuity, momentum, energy 
and magnetic hux. In the anelastic approximation (Gough 1969; Braginsky & Roberts 
1995; Lantz & Fan 1999; Jones et ah 2011) these equations take the form 


V • pn = 0, V ■ B = 0, 

R 5 1 

dfU + (y X u) X u = —Vn — r(fc X n) + —^r + F^ + -(V x B) x B, 


dtS + u-VS = 


P 

1 — Cl Pr / 1 

—V • RpTVS + ^{q, + 


FipT ' ' RT V"" ' Pmp 

dtB = V X {u X B) + Pm“^V^S, 


(2a) 

(2b) 

(2c) 

(2d) 


where u is the velocity, B is the magnetic hux density, S is the entropy and VII includes 
all terms that can be written as gradients. The viscous force (-Fj/), and the viscous {Qu) 
and Joule {Qj) heating are dehned in terms of the deviatoric stress tensor (Sij) 


F, = ^v-S, Q, = S:e, Q, = (V x Bf, (3) 

P 

where the double-dots symbol (:) denotes the Frobenius inner product. We assume that the 
viscosity z/ and the entropy dihusivity k are constant throughout the shell. The governing 
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equations are parametrized using the thickness of the shell d = Vo — Vi as a unit of length, 
d? jv as a unit of time, AS* as a unit of entropy, v^/jL^/d as a unit of magnetic induction, 
pc as a unit of density and Tc as a unit of temperature. Here, r* and Tq are the inner 
and the outer radius, A and po are the magnetic diffusivity and permeability, respectively. 
The system is then characterized by eight non-dimensional parameters: the radius ratio 
7] = ri/vo, the polytropic index n, the density scale number Np = In (p(rj)/p(ro)) , the 
Rayleigh number R = ciTcd‘^AS/{uK), the thermal Prandtl number Pr = z//k, the magnetic 
Prandtl number Pm = v/\, and the Coriolis number r = 2VLd?/v. 

Since the mass flux v = pu, and the magnetic flux density B are solenoidal vector 
helds, we employ a decomposition in poloidal and toroidal components. 


pn = V X (V X rrv) + V x rr'^w, 
B = 'V X (V X rh) -I- V X rg, 


(4a) 

(4b) 


where r is the radial unit vector, r is the length of the position vector r, v, w, h and g 
are the poloidal and toroidal scalars of the momentum and magnetic held, respectively. 
Equations (2a) are then satished automatically. Scalar equations for v and w are obtained, 
and effective pressure gradients are eliminated by taking f ■ V x Vx and r • Vx of 
Equation (2b). Similarly, equations for h and g are obtained by taking f ■ Vx and r- 
of Equation (2d). Spectral projections of the resulting poloidal-toroidal equations are 
presented in Appendix A, see also Section 2.3. 


2.2. Boundary conditions 


Equations (2) must be supplemented by boundary conditions. For the simulation 
results presented in this work the following boundary conditions are used. The inner and 
the outer surfaces of the shell are assumed as stress-free, impenetrable boundaries for the 
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flow 


P P 

f = 0, d^v ——dr{rv) = Q, drW —= 0, at r = Vi^Vo. (5) 

pr p 

In several cases discussed in Section 4 we use no-slip condition at the inner boundary, 

n = 0, drV = 0, tc = 0, at r = rj. (6) 

A fixed contrast of the entropy is imposed between the inner and the outer surface 

S' = 1 at r = Tj, S' = 0 at r = Tq. (7) 


The boundary conditions for the magnetic held are derived from the assumption of an 
electrically insulating external regions. The poloidal function h is then matched to a 
function which describes an external potential held, 

g = 0, h — = 0, dr{h — = 0 at r = r^, Tq. (8) 

We remark that our code presented below allows various other choices of boundary 
conditions to be made for all of the dynamical variables. 


2.3. Numerical method 

To perform the numerical simulations of this study, we have extended our Boussinesq 
code (Tilgner & Busse 1997; Busse et al. 2003; Simitev & Busse 2005; Busse & Simitev 
2006, 2008; Simitev & Busse 2009, 2012a) to solve the anelastic equations described in 
Section 2. Despite similarities with the Boussinesq code, this is a major modihcation both 
in terms of the mathematical model and the numerical code. For the numerical solution of 
the problem we have adapted the pseudo-spectral method described by Tilgner (1999). The 
scalar unknowns n, tc, h, g and S', are expanded in Chebychev polynomials Tp in the radial 
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direction r, and in spherical harmonics in the angnlar directions {0,(p) e.g., 

Nl I Nr 

V = ^ricosO) exp{imcp), (9) 

1=1 m=—l p=0 

where Pl^ denotes the associated Legendre fnnctions, x{r) = 2(r — r*) — 1, and Ni and 
are trnncation parameters. A system of eqnations for the coefficients in these expansions is 
obtained by a combination of a Galerkin spectral projection of the governing eqnations in 
the angnlar directions and a collocation constraint in radius. This system is presented in 
Appendix A. 

Computation of nonlinear terms in spectral space is expensive, so nonlinear products 
and the Coriolis term are computed in the physical space and then projected into the 
spectral space at every time step. A standard 3/2-deahasing in 9 and (p is used at this 
stage. A hybrid of a Crank-Nicolson scheme for the diffusion terms and a second order 
Adams-Bashforth scheme for the nonlinear terms is used for integration in time. 

Calculations are considered adequately resolved when the spectral power of kinetic 
and magnetic energy drops by more than two orders of magnitude from the spectral 
maximum to the cut-off wavelength as suggested e.g. by Christensen et ah (1999). A range 
of numerical resolutions has been employed in this study varying from {Nr = 61, A') = 96) 
in less demanding cases to {Nr = 121, A) = 144) in more strongly stratihed or turbulent 
runs. Correspondingly, the physical gridpoints on which non-linear terms are evaluated 
have been varied up to Nr = 121, No = 216, N^p = 437. We hnd that this provides adequate 
resolution as demonstrated in Figure 1 for a typical dynamo solution. 

The pseudo-spectral approach described above is the most common method for solving 
the fundamental equations of convection-driven flow and electromagnetic induction within 
a rotating spherical shell hlled with an electrically conducting fluid. The approach was 
pioneered by Clatzmaier (1984) and with appropriate modihcations it has been widely used 
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by various groups for modelling convection-driven geo-, Solar, and planetary dynamos. A 
number of codes based on similar principles have been developed that differ mainly in details 
such as time stepping methods and treatment of radial dependence with hnite differencing 
and Chebyshev decomposition in r being two popular choices. Early versions of some codes 
were derived directly from the code of Glatzmaier (1984), including the anelastic ASH 
code extensively used for Solar simulations (Clune et al. 1999) and the Boussinesq MAG 
code used for geodynamo simulations (Olson et al. 1999). Other codes including ours were 
developed independently (e.g. Jones et al. 1995; Tilgner 1999; Hollerbach 2000). While it is 
not feasible to provide here a comprehensive list of existing numerical codes and discuss the 
numerous variations in actual implementation, we refer to a series of benchmarking papers 
(Christensen et al. 2001; Jones et al. 2011; Marti et al. 2014; Jackson et al. 2014) and to the 
reviews (Miesch 2005; Wicht & Tilgner 2010) for overview of codes commonly used in the 
solar context and in the geo-/planetary context, respectively. These references also include 
discussions of other essentially different numerical approaches of solution. Recent studies 
related to ours (Fan & Fang 2014; Karak et al. 2015; Mabuchi et al. 2015; Hotta et al. 2015) 
use hnite difference methods to solve comparable but different sets of equations. Local 
methods have better parallel efhciency but also inferior accuracy (Tilgner 1999). They also 
have difhculties with imposing global boundary conditions for the magnetic held, e.g. all of 
the above dynamo models use the unphysical radial condition for the magnetic held on the 
outer surface. They also encounter difhculties in treating spherical geometries, e.g. Fan & 
Fang (2014) and Karak et al. (2015) consider wedges and not spherical shells. Our code has 
been developed and optimized independently over a number of years and a large database 
of Boussinesq results is available for comparison. The current anelastic version of our code 
is new and perhaps unique among other spectral codes, with the exception of the ASH 
code, in allowing for a radial dependence of the viscosity and the thermal and magnetic 
dihusivities. However, the latter facilities have not been used in the present analysis. 
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2.4. Diagnostic output quantities 


We characterise convection and dynamo solutions by their kinetic and magnetic energy 
and heat transport given by a Nusselt number. The energies can be conveniently split 
into poloidal and toroidal components, mean and fluctuating components and further into 
equatorially-symmetric and equatorially-antisymmetric components. We thus obtain a good 
characterisation of the scales of the convective flow and the multipole structure of dynamos. 
The mean and fluctuating toroidal and poloidal components of the total kinetic energy E'kin 
are defined as 


Ep = {{Vx{Vvxr)Y/{2p)), E, = {{Vriu x r)"/{2p)), (10a) 

Ep = ((V X (Vh X r))V(2p)), Et = {{Vrw x rf/{2p)), (10b) 


where angular brackets {) denote averages over the spherical volume of the shell, overlaid 
lines denote axisymmetric parts and overlaid check marks denote non-axisymmetric parts 
of a scalar held. The total magnetic energy T^magn can be split in a similar way with 
components defined as in Equations (10) but with h and g replacing v and w and without 
the factor p~^ within the angular brackets. The total energies are, of course, the sum of 
all components. The Nusselt number is defined as the ratio between the values of the 
luminosity of the convective state Lconv and the luminosity of the basic conduction state 

L basic; 


Nu = 


L 


Lprinv 


'basic 


fdV{r) 


KpT{drS)r'^ sin 9d9dip, Lbasic = 


47rncip(r)"' 

exp(At^)-l’ 


with the integral taken over a spherical surface dV (r) at radius r. Apart from quantifying 
the heat transport of convection, the value of the Nusselt number serves as a convenient 
proxy for the super-criticality of the convective regime. 


Other diagnostic quantities used below are the non-dimensional magnetic Reynolds 
R. . Rp... ..P., Ro . 
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2.5. Benchmarking and validation 

To validate the new code we present a comparison with the anelastic benchmark 
simulations recently proposed by Jones et ah (2011). For the comparison we employ an 
alternative parametrization based on the magnetic rather than viscous diffusion time scale, 
used in the benchmark models. Our output results from the three benchmark cases defined 
in Jones et ah (2011) are summarized in Table 1, and selected components of the solution 
are plotted in Figure 2. The mean values and the means of the deviations computed from 
the values reported by the four codes participating in (Jones et ah 2011) are also listed in 
Table 1. We achieve perfect agreement with the benchmark results for the hydrodynamic 
case and the steady dynamo case, labelled Benchmark 1 and Benchmark 2. Our results for 
the unsteady dynamo case labeled Benchmark 3 show some insignificant differences from 
the values reported in Jones et al. (2011). The reason for the discrepancies is the use of 
imposed 2-fold azimuthal symmetry and lower resolution in our code to reduce computing 
time. 


3. Differential rotation and dynamo action in the buoyancy-dominated regime 

It is well known, that as buoyancy forcing becomes significantly larger than the Coriolis 
force anti-solar differential rotation develops. Such buoyancy-dominated regime was first 
identified by Gilman (1976) and Gilman & Foukal (1979), and more recently it was studied 
by Aurnou et al. (2007) and Gastine et al. (2013). More recent studies include the works of 
(Guerrero et al. 2013a,b) that is closely tailored to the solar case. These studies consistently 
found that due to vigorous mixing angular momentum is homogenized within the whole 



volume of the shell, and this leads to a strong retrograde zonal flow in the equatorial region, 
and thus to the anti-solar type of rotation prohle. 


3.1. Transition between rotation-dominated and buoyancy-dominated regimes 

Here, we investigate the transition from the rotation-dominated regime to the 
buoyancy-dominated regime for non-magnetic and magnetic (dynamo) cases. Figure 3 
demonstrates the transition from the solar-type to the anti-solar type of differential rotation 
in a set of cases with increasing value of the Rayleigh number R and all other parameter 
values kept hxed. The increase of the Rayleigh number is the most direct approach to 
the buoyancy-dominated regime as R is a direct measure of the magnitude of buoyancy 
forcing. We remark, however, that the buoyancy-dominated regime can also be reached if 
other control parameters are varied. For instance, a decrease in the Coriolis number r at 
hxed values of the other parameters including R will equally well bring convection to the 
buoyancy-dominated regime since r is a measure of the Coriolis force. Reaching the regime 
by variation of the other parameters is also possible but less straightforward. In the case 
of Figure 3 the transition to the buoyancy-dominated regime happens between R = 10^ 
and R = 2 X 10^. The Coriolis-dominated regime is characterized by columnar convection 
structures that are oriented parallel to the rotation axis and mostly outside the tangent 
cylinder, the cylinder that touches the inner core at the equator. 

In the case of Figure 3, columnar convection has a dominant azimuthal wave number 
of 8. This number, however, varies as the parameter values are varied. Since even in the 
strongly chaotic regime the azimuthal wave number of convection remains similar to that 
near the onset of convection, it is useful to note the study of Busse & Simitev (2014) 
where the critical onset for this problem has been studied. Differential rotation in the 
Coriolis-dominated regime is solar-like and geostrophic, i.e. constant on cylindrical surfaces. 
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Meridional circulation is relatively weak but changes from single to two cells in a hemisphere 
with the increase of R. In this connection we note recent observational results by Zhao 
et al. (2013) that report two-cell meridional circulation in the solar convection zone. 

Convection in the buoyancy-dominated regime for R > 2 x 10^ in the case of Figure 
3 becomes disorganized. Convective columns are broken, and the convection pattern loses 
its anisotropy with respect to the axis of rotation. Convection, for instance, is now not 
restricted to the region outside of the tangent cylinder but produces vigorous flows in the 
polar regions as well. Meridional circulation appears to return to a single-cell pattern but 
the symmetry with respect to the equatorial plane is lost. This, of course, is due to the fact 
that the Coriolis force is no longer dominant and the role of rotation is much diminished. 

The most notable effect is the sign reversal of the differential rotation that switches 
from the solar-like to the anti-solar prohle. The anti-solar differential rotation remains 
largely constant on cylinders parallel to the rotation axis. While some conical features are 
present in the case of Figure 3 they seem to be conhned near the surface and appear less 
signihcant. 


3.2. Effects of magnetic field on differential rotation 

While solar convection is likely dominated by buoyancy, it is rather challenging to 
reconcile the anti-solar differential rotation found in the buoyancy-dominated regime with 
observations. It is well known from helioseismology observations that solar differential 
rotation is prograde and strongly non-geostrophic (e.g. Thompson et al. 1996). In this 
situation, magnetic effects may provide one possible mechanism for reversing the anti-solar 
differential rotation into the solar-like type (Fan & Fang 2014; Karak et al. 2015; Mabuchi 
et al. 2015). This is not unreasonable to expect. Indeed, it is well-established that the main 
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effect of self-generated magnetic field on convection is to impede differential rotation. This 
effect is measured by a strong decrease of mean toroidal kinetic energy in dynamo solutions 
initially demonstrated by Grote & Busse (2001) and later studied in the parameter space 
by Simitev & Busse (2005) (see their Fig. 16) and Busse & Simitev (2005) (see their Fig. 
7). It has also been investigated by Aubert (2005) and Yadav et ah (2013). To explore the 
effect of magnetic held on buoyancy-dominated convection as described in the preceding 
Section, we have performed a systematic comparison between large sets of simulations of 
non-magnetic convection and of self-sustained dynamos. 

Figures 4 and 5 summarize the results of four such sequences of cases with one pair 
at Prandtl number Pr = 1 and Coriolis number r = 2000 and the other pair at somewhat 
smaller values Pr = 0.5 and r = 300. One sequence in each pair is non-magnetic while the 
other sequence includes a self-generated magnetic held. The specihc choice of parameters is 
motivated by properties of solar convection as discussed below. Time-averaged components 
of the kinetic and magnetic energy densities are plotted as a function of the increasing 
Rayleigh number R for both the convection-only cases and the dynamo cases. The top row 
of Figure 4 shows the time-averaged value of the diherential rotation at the equator on the 
outer surface of the spherical shell, M(^(ro, 7r/2), and serves as an easily accessible indicator 
of the type of diherential rotation. The transition from the rotation-dominated regime 
(solar-like rotation) to the buoyancy-dominated regime (anti-solar rotation) is thus easily 
identihed by the sign change of u^piro, vr/2) (marked by a dash-dotted line in Figures 4 and 
5). Convection in the absence of magnetic held is characterized by an abrupt as opposed 
to a gradual increase of both diherential rotation and meridional circulation. The total 
magnetic energy of the self-sustained held is on average an order of magnitude smaller than 
the kinetic energy, but the magnetic held has a signihcant ehect on both the diherential 
rotation and meridional circulation. In contrast to the non-magnetic case, they are strongly 
reduced and show no abrupt change in their values. 
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In the sequences shown in Figure 4(a,c,e) the transition from the solar-like to anti-solar 
differential rotation is, in fact, suppressed. While we expect that if buoyancy is further 
increased, i.e. by increasing R, convection will eventually arrive once again at the transition 
to prograde rotation, we emphasize that the suppression of the transition happens over 
a large interval of R comparable to the interval between the onset of convection and the 
solar-antisolar transition itself. 

In the sequences shown in Figure 4(b,d,f) the suppression of the solar-antisolar 
transition is not observed, even though the effects of decreasing differential rotation and 
meridional circulation are visible. This difference illustrates other important parameter 
dependences, primarily those on the Prandtl and the Coriolis numbers, Pr and r, 
respectively. In the region of low values of Prandtl and Coriolis numbers convection is 
known to be rather different from columnar convection in that it takes the form of an 
equatorial belt of large cells attached near the outer surface of the spherical shell (Ardes 
et al. 1997; Busse & Simitev 2004). Differential rotation generated by equatorially-attached 
convection is typically less affected by the braking effect of the magnetic held (Simitev & 
Busse 2005). 

Figure 5 shows the values of the Rossby, Lorentz and magnetic Reynolds number in 
dependence of the Rayleigh number of the same sequences as shown in Figure 4. The 
transition from rotation-dominated to buoyancy-dominated convection happens at about 
Ro = 1 in the sequence of Figure 4(a,c,e) and at about Ro = 0.5 in the sequence of Figure 
4(b,d,f). These values are similar to the ones reported by Gastine et al. (2014), Karak 
et al. (2015) and Mabuchi et al. (2015) but a weak dependence on the Prandtl and the 
Coriolis numbers appears to exist. The dynamo effects do not appear to affect the value of 
the Rossby number for transition in agreement with Karak et al. (2015) and Mabuchi et al. 
(2015). 
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Further notable effects of the magnetic field are summarized in Table 2. These effects 
are illustrated in terms of one selected strongly chaotic case discussed below. 


3.3. Structure and dynamics of convective flows and magnetic fields in the 

bnoyancy-dominated regime 

Figure 6 shows a comparison of the spatial structures of the convective flow of 
a non-magnetic case and those of a self-sustained magnetic dynamo case at identical 
parameter values 7] = 0.65, Pr = 1, r = 2000, R = 10^, n = 2, Np = 3 and for the dynamo 
solution Pm = 2. At rj = 0.65 the shell thickness is slightly thicker than the thickness 
of the solar convection zone and is selected for ease of numerical simulation. The typical 
size of convective structures is related to the thickness of the shell and thus thinner shells 
require spherical harmonics decomposition of higher order and degree to resolve the angular 
structure of the flow. The value Pr = 1 is appropriate in the sense that turbulent mixing 
tends to homogenize the flow and molecular diffusivities are replaced by effective turbulent 
diffusivities of similar magnitude. At r = 2000 the Coriolis number is moderately but 
not excessively large reflecting the model assumption that the flow in the deep convection 
zone is buoyancy rather than rotation dominated. The values of the polytropic index 
n = 2 is adequate, while the value of the density scale height Np = 3 is much lower than 
estimated for the solar convection zone. However, increasing, Np much beyond 5 becomes 
computationally very demanding. Finally, the value of the Rayleigh number has been 
selected such that the non-magnetic convection case is located in the buoyancy-dominated 
regime, the onset of which is at R = 5.2 x 10® for these parameter values. 

In the case of non-magnetic convection shown in the right column of Figure 6 differential 
rotation is in anti-solar direction. It is strongly geostrophic i.e. constant on cylinders 
parallel to the rotation axis. Differential rotation is monotonously increasing towards the 
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outer surface of the spherical shell. Note that this is also true for solar-like differential 
rotation in the rotation dominated regime. The structure of the flow changes signihcantly 
with radius. The hrst two rows of Figure 6 show isocontours of the radial velocity near 
the surface and somewhat below mid depth within the spherical shell. The flow near the 
surface is a patchwork of small-scaled up- and down-wellings distributed in a very chaotic 
pattern over the full surface of the spherical shell. No visible structure can be discerned 
and the location of the convective cells changes chaotically in time (not shown). The scale 
of the convective structures increases in depth and near the equator elongated convective 
cells tilted clockwise in the northern hemisphere and tilted anticlockwise in the southern 
hemisphere are found. On average this equatorial pattern drifts in the retrograde direction 
carried away by the strong anti-solar differential rotation. The so described radial structure 
is also evident in the contour lines of the radial velocity in the equatorial plane plotted in 
the third row of 6. Finally, the meridional flow takes the form of two large circulations 
in poleward direction at the surface of the shell. The meridional circulations are nearly 
strictly mirror-symmetric with respect to the equatorial plane. This and the symmetry of 
the differential rotation are remarkable large-scale coherent features of this otherwise very 
chaotic solution. 

As noted above the magnetic energy is signihcantly lower than the kinetic energy of the 
how. The inhuence of the magnetic held on convection, however, in the dynamo case shown 
in the left column of Figure 6 is quite remarkable. The most notable ehect is, of course 
the reversal of the direction of diherential rotation from the anti-solar to solar-like type. A 
further remarkable diherence is that the maximum of the diherential rotation occurs in the 
depth of the spherical shell rather than at the surface as is always observed in the case of 
non-magnetic convection. This is potentially a signihcant ehect as it means that there is a 
negative gradient of diherential rotation in the subsurface layer of the shell. 
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In the magnetic case the structure of the flow also changes signihcantly with radius. 
The hrst two rows of Figure 6 show isocontours of the radial velocity at the same radial 
values as in the non-magnetic case. The flow near the surface is a again a patchwork 
of small-scaled up- and down-wellings distributed in a very chaotic pattern. No visible 
structure can be discerned and the location of the convective cells changes chaotically 
in time (not shown). An important effect of the magnetic held is that in the magnetic 
dynamo case convection appears to be stronger in the polar regions rather than outside the 
tangent cylinder. The main difference, however, appears in depth. Large scale convective 
columns arranged in a cartridge belt pattern within the tangent cylinder and spanning both 
hemispheres to about mid latitudes are clearly visible. The columns drift in the prograde 
direction due to the solar-like differential rotation. This is also a signihcant observation 
because it indicates that very little may be inferred for the structure of deep convection 
from observations of near surface hows. In particular, it is not known whether large scale 
convective columns exist in the deep solar convection zone or not. Our results indicate 
large scale convective columns hidden from view by much smaller-scale chaotic convection 
with no discernible structure is a likely dynamical possibility. Finally, the meridional how 
of the magnetic dynamo case appears rather disorganized with a number of smaller scale 
circulations appearing in both hemispheres as shown in the bottom row of Figure 6. 

The structure of the generated magnetic held in the dynamo case is shown in Figure 
7. The magnetic held has a large-scale dipole component emerging from a patchwork of 
small scale magnetic features. The dipole is mainly supported by strong polar magnetic hux 
tubes in the polar region, which in turn are due to the relatively strong polar convection. 
The predominant polarity is less clear in the equatorial region where the magnetic held 
structures are smaller in scale and of both polarities. The dipole solution is non-oscillating. 
Unfortunately, we have not been able to locate oscillating dynamos in this regime and to 
observe the direction of dynamo wave propagation. This is left for future studies. 
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4. Remarks on anelastic dynamos in the rotation-dominated regime 

While the attention in this paper is focussed on dynamo effects near the transition 
from rotation-dominated to buoyancy-dominated convection, in this section we wish to 
demonstrate that conical prohles of differential rotation as well as regular and persistent 
dynamo oscillations can be found in the parameter space of our minimal convection-driven 
dynamo model without recourse to additional modelling assumptions. We also take the 
opportunity to elucidate some of the points made about buoyancy-dominated dynamos 
above by comparison with features in the rotation-dominated regime. 

Perhaps the only clear example of a conical prohle of differential rotation in a single¬ 
layer simulation with spherically-symmetric boundary conditions is that reported in the 
work of Brun & Toomre (2002). However, these authors use subgridscale parametrisation 
of diffusivities which clearly affects results. Following Miesch et ah (2006), the majority of 
models that report conical prohles seem to impose a non-zero latitudinal gradient of entropy 
as their bottom boundary condition (e.g. Fan & Fang 2014) or to include anisotropic 
heat conductivity (e.g. Karak et al. 2015) or a stably stratihed layer at the bottom of the 
convection zone (e.g. Mabuchi et al. 2015) all of which increase the baroclinicity and induce 
conical prohles. In this context. Figure 8 shows an example of diherential rotation with 
some conical features in the lower part of the convection zone, obtained in our minimal 
self-consistent formulation of the problem. The parameter values of this run are the same 
as those for the sequence of cases reported in Figure 4(a,c,e), with a value of the Rayleigh 
number that places it in the rotation-dominated regime, and a somewhat larger value of 
Pm which is known to promode stronger dipolar helds (Simitev & Busse 2005). This run 
dihers from the latter sequence only in that it uses a no-slip velocity condition on the inner 
spherical boundary. Convection in the polar regions is weak if not fully absent. Convective 
hows are conhned outside of the tangent cylinder and this is where the dynamo process is 
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also located resulting in a magnetic field that is strong near the equator and at midlatitudes 
but weak in the polar regions. We wish to contrast this situation with the situation 
discussed in connection to the dynamo case presented in Figures 6 and 7. This comparison 
makes it rather obvious that in the latter case vigorous convection in the polar regions 
gives rise to strong magnetic held in the same regions. The dynamo shown in Figure 8 is 
an oscillatory dynamo and the comparison with the case of Figure 7 elucidates the reasons 
why the latter is non-oscillatory. While in the case of Figure 8 both the magnetic held and 
the diherential rotation achieve their maxima in the same region (the equatorial region), in 
the case of Figure 7 the regions where the maximal amplitude of the magnetic held and 
of the diherential rotation occur do not coincide which is detrimental to aVL oscillations 
(Busse & Simitev 2006; Warnecke et ah 2014). 

To illustrate the oscillations in question, we present in Figure 9 one period of a 
predominantly dipolar dynamo wave. The parameter values of this run are identical to 
the cases shown in Figure 8 except for a slightly larger value of the Rayleigh number 
which helps to make the oscillations more regular. The dynamo wave is driven by the 
aVL mechanism hrst proposed by Parker (1955) and later conhrmed in three-dimensional 
simulations by Busse & Simitev (2006), see also (Simitev & Busse 2012a; Schrinner et al. 
2012; Warnecke et al. 2014). The dynamo wave propagates in the direction of the poles. 
This case shows similar conical features of the diherential rotation prohle, and Figure 9 also 
illustrates their variations in time. Figure 10 shows the dominant dipolar and quadrupolar 
components of the magnetic held represented by time series of the appropriate coefficients 
in the spherical harmonic expansions of the toroidal and poloidal scalars of the magnetic 
held. The time series show that the oscillations are very regular and persistent over the 
course of the simulation and that the dipole is dominant. 

We wish to conclude this section with the remark that the parameter space of the 


minimal model formulation seems to merit further investigation especially in the case when 
self-generated magnetic helds are present. 


5. Conclusion 

We have presented in this paper several sets of convective dynamo simulations in 
density-stratihed rotating spherical fluid shells based on a physically consistent anelastic 
model with a minimum number of parameters. The computations are performed using 
a new simulation code that is also presented here for the hrst time along with code 
validation results against published benchmark solutions. We demonstrate that conical 
differential rotation prohles and persistent regular dynamo oscillations can be obtained in 
the parameter space of this minimal formulation without recourse to additional modelling 
assumptions. The main focus of the work is placed on extending the dynamo simulations 
into a “buoyancy-dominated” regime where the buoyancy forcing is dominant while the 
Coriolis force is no longer balanced by pressure gradients and where strong anti-solar 
differential rotation develops as a result. The dynamo solutions are compared to identical 
sets of non-magnetic convection solutions to reveal the effects of the self-sustained magnetic 
held on convection in general and on the differential rotation in particular. The most 
signihcant results are summarized in Table 2 and below we discuss some similarities between 
our solutions and solar and stellar convection. 

We also wish to compare our results to studies on a similar topic reported in the 
recent literature, in particular the works of Fan & Fang (2014), Karak et al. (2015) and 
Mabuchi et al. (2015) where dynamo simulations are reported and of Gastine et al. (2013) 
and Hotta et al. (2015) where hydrodynamic simulations are reported. Before we comment 
on similarities and differences in solutions, we wish to point out that with the exception 
of the model considered by Gastine et al. (2013), the models considered by the other 
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groups have significant and essential differences compared to ours. The models of Hotta 
et ah (2015), Karak et ah (2015) and Mabuchi et ah (2015) are fully compressible models 
in contrast to our anelastic approximation. Hotta et ah (2015) use artihcially enhanced 
viscosity, a radius dependent cooling term and are interested mainly in the properties of the 
near-surface convection layer. The model of Mabuchi et ah (2015) consists of two layers - 
a stably stratihed layer surrounded by a convective envelope. The models of Fan & Fang 
(2014) and Karak et al. (2015) consider wedges, i.e. partial spherical shells, and are not 
fully spherical, so polar convection is effectively not represented in their solutions. This 
partly explains why they hnd oscillatory dynamos. Conical differential rotation prohles 
are promoted in the latter four models due to the inclusion of secondary physical effects: 
Fan & Fang (2014) impose a latitudinal entropy gradient as a bottom boundary condition. 
Fan & Fang (2014), Hotta et al. (2015) and Karak et al. (2015) include radial variation 
in diffusivities and parametrisations of unresolved scales. Numerical implementations also 
differ - all codes except that of Gastine et al. (2013) are based on hnite-difference methods 
while ours is pseudo-spectral. An artihcial radial magnetic held condition is imposed on the 
outer boundary in these studies which may signihcantly distort dynamo effects. Despite 
the differences in modelling strategy, it is signihcant that we hnd a number of similarities 
in our simulations. This increases the conhdence in the robustness of the reported results. 

Oscillations obtained in dynamos generated within the rotation-dominated regime with 
few exceptions (Warnecke et al. 2014), appear to always travel in the poleward direction 
much like illustrated in Figure 9. For references see (Busse & Simitev 2006; Simitev & 
Busse 2012a), and also (Schrinner et al. 2012). There is some evidence that dynamo waves 
can travel towards the equator when a negative radial gradient of the diherential rotation 
prohle exists (Simitev & Busse 2012b), and when in addition to the latter the a-ehect, 
proportional to —(V x u) ■ u + {V x B) ■ B/p, is positive (negative) in the northern 
(southern) hemisphere (Warnecke et al. 2014). This evidence supports an early analysis by 
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Yoshimura (1975). While we commonly find dipolar oscillations in the rotation-dominated 
regime, our dynamos in the buoyancy-dominated regime do not oscillate. This is in contrast 
to the results of Fan & Fang (2014) and Karak et ah (2015) which may be attributed to the 
absence of polar convection in their simulations due to the absence of conical polar sections 
used in their geometrical configuration. On the other hand, we agree with Mabuchi et ah 
(2015) who use a full spherical shell and find that dynamos with anti-solar rotation are 
predominantly dipolar and non-oscillatory. Despite being relatively weak the self-sustained 
magnetic fields in the buoyancy-dominated regime reported in the present study are able 
to reverse the direction of differential rotation from anti-solar to solar-like. From the 
perspective of oscillations, it is significant that we find that differential rotation attains a 
maximum inside the shell and that a negative radial gradient is persistently maintained in 
the near-surface layer. This may facilitate equatorward dynamo wave propagation in the 
buoyancy-dominated regime. The differential rotation in our buoyancy-dominated dynamos, 
e.g. Figure 6, has a cylindrical profile, in contrast to the more conical profiles reported by 
Fan & Fang (2014) Karak et ah (2015) and Mabuchi et ah (2015). This difference is almost 
certainly caused by the fact that a non-zero latitudinal gradient of entropy is imposed as a 
bottom boundary condition by Fan & Fang (2014), that an anisotropic heat conductivity 
is used by Karak et al. (2015), and that a stably stratified layer at the bottom of the 
convection zone is present in (Mabuchi et al. 2015). 

We find that the dynamo-generated magnetic field can suppress the transition from 
the solar-like to the antisolar-like rotation profile thus confirming similar findings reported 
by Fan & Fang (2014) Karak et al. (2015) and Mabuchi et al. (2015). In this case the 
convection is significantly stronger near the poles than in the equatorial region, leading 
to a predominantly dipolar dynamo where both the toroidal and the poloidal fluxes are 
stronger in the polar regions compared to equatorial regions. Such dynamo regime with 
concentration of magnetic field in the polar regions may explain the observations of polar 
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starspots in young solar-type stars which exhibit reduced but still the solar-type differential 
rotation (Brown et ah 2014), also see (Yadav et ah 2015). While starspots reduce locally 
the vigour of convection in the near surface layer of the convection zone, sufficiently strong 
convection in depth is required to generate magnetic helds that are large enough to cause 
starspots in the hrst place. Our calculations do not have sufficient resolution to resolve 
strongly turbulent stellar near surface layers. 

Our simulations conhrm the hndings of Gastine et al. (2013) and Hotta et ah (2015) 
that different regimes of convection occur in the inner and at the outer part of the 
spherical shell simultaneously such that organized geostrophic convection columns are 
hidden below a near-surface layer of well-mixed highly-chaotic convection. Both of the 
latter studies are non-magnetic and the work of Hotta et al. (2015) reports simulations in 
the rotation-dominated regime only. The model of Gastine et al. (2013) is rather similar 
to ours and it is one of the aims of the present paper to extend their analysis through 
considerations of dynamo effects. On the Sun small scale turbulent convection is clearly 
observable in the subsurface layer of the solar convection zone while simulations inevitably 
hnd some columnar structures. Evidence of different convection morphology as a function 
of radius is significant because it provides a bridge between observations and simulations. 
The deeper large-scale organized convection columns are likely to play important role in the 
solar dynamo and its magnetic cycles. 
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A. Spectral projection of the toroidal-poloidal governing equations 


Scalar equations for v and w are obtained, and effective pressure gradients are 
eliminated by taking f ■ V x Vx and f ■ Vx of Equation (2b). Similarly, equations for 
h and g are obtained by taking f ■ Vx and r- of Equation (2d). The scalar unknowns 
V, w, h, g and S are then expanded in Chebychev polynomials Tp in the radial direction 
r, and in spherical harmonics 1)™ in the angular directions {6, ip) as shown in Equation 
(9). After a standard Galerkin projection procedure in the angular directions {0,(p) 
the following set of partial differential equations for the spectral expansion coefficients 
{u["(r, t),w‘^{r, t), h‘^{r, t),g'j^{r, t), S]^{r, t) : 1 = 1 ... m = —/.../} is obtained 
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and the following notation of some radial functions is used for brevity 
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The expression for the viscous dissipation in Equation (Ale) is 


Q, = Qa)+Q(2), 
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Finally, square brackets with a subscript (and superscript) as in or [-j™ denote a component of 
a vector or an appropriate coefficient in a spherical harmonic expansion, respectively. 


Solution of Equations (Al) proceeds as described in Section 2.3. 
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Fig. 1.— Time-averaged power spectra of kinetic (full circles) and magnetic (empty circles) 
energy as a function of (a) the harmonic degree I, (b) the harmonic order m and (c) the 
Chebyshev polynomial degree p in the case p = 0.65, R = 6 x 10®, Pr = 1, Pm = 2, r = 2000, 
n = 2, Np = 3. 
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Table 1: Comparison with the benchmark solutions reported by Jones et ah (2011). The 
values labelled “Mean ± Mean Deviation” are the respective means and the means of the 


deviations computed from the values reported by the four codes participating in (Jones et ah 
2011 ). 



Benchmark 1: 

Hydrodynamic 

convection 

Benchmark 2: 

Steady dynamo 

Benchmark 3: 

Unsteady dynamo 

h 

0.35 

0.35 

0.35 

n 

2 

2 

2 

Np 

5 

3 

3 

Pr 

1 

1 

2 

Pm 

1 

50 

2 

r 

2000 

1000 

4 X 10^ 

R 

351806 

8 X 10^ 

2.5 X lO'^ 

Nr / Nr 

129 / 129 

129 / 129 

111 / 111 

Ni / Ne 

128 / 128 

128 / 128 

120 / 144 

Nm/ % 

129 / 257 

129 / 257 

61 / 73 

Timestep 

4 X 10"® 

1 X 10“® 

1 X 10"'^ 

E 

81.87991 

4.19405 X 10® 

2.32730 X 10® 

Mean ± Mean Dev. 

81.680 ± 0.245 

(4.186 ±0.013) X 10® 

(2.317 ±0.014) X 10® 

Ep 

0.02201 

53.0100 

100.40 

Mean ± Mean Dev. 

0.0220 ± 0.0001 

52.90 ±0.15 

111.75 ± 3.75 

Et 

9.37598 

6.01725 X 10^ 

1.81399 X 10"^ 

Mean ± Mean Dev. 

9.3568 ± 0.0282 

(6.001 ±0.018) X 10"^ 

(1.355 ±0.008) X 10^ 

M 

Mean ± Mean Dev. 

- 

3.20172 X 10® 

(3.194 ±0.088) X 10® 

2.58012 X 10® 

(2.413 ±0.023) X 10® 

Mp 

Mean ± Mean Dev. 

— 

1.69650 X 10^ 

(1.692 ±0.005) X 10"^ 

2.91155 X 10^ 

(2.155 ±0.070) X 10^ 

Mt 

Mean ± Mean Dev. 


2.41185 X 10® 

(2.412 ±0.028) X 10® 

1.17292 X 10^ 

(0.948 ± 0.003) X 10^ 

Luminosity 

4.19886 

11.50302 

42.50992 

Mean ± Mean Dev. 

4.19886 ±3 X 10-® 

11.503 ±4 X 10“® 

42.75 ±0.15 
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Table 2: Summary of the effects of self-sustained magnetic field and comparison with non¬ 
magnetic convection. 


Non-magnetic convection 

Dynamo 

Monotonic increase/decrease of differ- 

Differential rotation attains a maximum 

ential rotation towards the outer surface 

inside shell and a subsurface decrease. 


in solar/antisolar cases. 


Retrograde differential rotation in 

Differential rotation reversed from anti- 

buoyancy-dominated regime. 

solar to solar-like. 

No columnar structure at depth. 

Convective columns visible in depth. 
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Benchmark 1: Benchmark 2: Benchmark 3: 



Fig. 2.— Solutions structures of benchmark cases 1, 2 and 3 (left to right). The first plot in 
each column shows azimuthally-averaged isocontours oiu^p (left half) and of the streamlines 
r sin 6^(501;) (right half) in the meridional plane. The second plot in each column shows 
isocontours of Ur in the equatorial plane. The third plot in each column shows isocontours 
of Br aX r = To- 
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Fig. 3.— Structures of convection showing the transition to the buoyancy-dominated regime 
with increasing value of the Rayleigh number as indicated in the plot and rj = 0.65, Pr = 0.3, 
Pm = 3, r = 200, n = 2 and Np = 3. The plots in the hrst column show time- and 
azimuthally-averaged isocontours oiuip (left half) and of the streamlines r sin 6*((96iu) (right 
half) in the meridional plane. The plots in the second column show contours of instantaneous 
Ur on the spherical surface r = {ri + r’o)/2. The plots in the third column show contours of 
instantaneous Ur in the equatorial plane. 
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Fig. 4.— (a,b) Differential rotation at the equator 7r/2). (c,d) Average kinetic energy 

densities and (e,f) average magnetic energy densities as functions of the Rayleigh number R 
in the cases (a,c,e) for rj = 0.65, Pr = 1, r = 2000, n = 2, Np = 3, and in the cases (b,d,f) 
for T] = 0.65, Pr = 0.5, r = 300, n = 2, Np = 3. Nonmagnetic convection cases are denoted 
by thin symbols in (a,b,c,d). Dynamo cases are denoted by thick symbols in all panels and 
have Pm = 2 in (a,c,e) and Pm = 6 in (b,d,f). Black circles, red squares, green pluses and 
blue crosses denote Xp, Xt, Xp, X^, with X = E, M. Vertical dash-dotted lines denote the 
transition to buoyancy-dominated regime. 
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Fig. 5.— Time-averaged values of the Rossby number Ro (black diamonds) and the Lorentz 
number Lo (red circles) measured on the left-hand axis and of the magnetic Raynolds number 
Rm (blue crosses) measured on the right-hand axis. Panel (a) shows the same sequences 
illustrated in Figure 4(a,c,d) and panel (b) shows the same sequences illustrated in Figure 
4(b,d,e). The thin symbols indicate the values of Ro for non-magnetic convection and bold 
symbols represent the dynamo cases. 
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Fig. 6.— Comparison between non-magnetic convection (left) and dynamo (right) at iden¬ 
tical parameter values r] = 0.65, Pr = 1, r = 2000, R = 10^, n = 2, Np = 3 and for the 
dynamo Pm = 2. Plots in the hrst row show Ur at r = 0.95 -|- r*, plots in the second row 
show Ur at r = 0.3 -|- rj, plots in the third row show Ur in the equatorial plane, plots in 
the fourth row show isocontours of the differential rotation (left half) and of the stream¬ 
lines r sm6{dgv) (right half) in the meridional plane. The plots in the hrst three rows are 
instantaneous snapshots, while the density plots in the fourth row are time-averaged. 
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Fig. 7.— Magnetic field components of the same case as in the right column of Figure 6. 
The left plot shows contours of Br at r = 1.13 + rj, the right plot shows contours and 
meridional held lines. 
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Fig. 8.— Flow and field structures in the case 77 = 0.65, Pr = 1, r = 2 x 10^, R = 1.8 x 10®, 
Pm = 8, n = 2, Np = 3 with no-slip condition on the inner boundary. The top left 
plot shows time-averaged isocontours of the differential rotation (left half) and of the 
streamlines r sin 9{dev) (right half) in the meridional plane. The top right plot shows contours 
of instantaneous at r = 0.5 -|- r^. The bottom left plot shows contours of instantaneous 
B^p (left half) and meridional field lines (right half). The bottom right plot shows contours 
of instantaneous Bj. at r = 1.13 -|- r*. 
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Fig. 9.— One period of dipolar oscillations in the case 77 = 0.65, Pr = 1, r = 2 x 10^, 
R = 2.5 X 10®, Pm = 4, n = 2, iVp = 3 with no-slip condition on the inner boundary. Time 
between plots is At = 0.09 staring at to = 96.8174 in the time series shown in hgure 10. 
The hrst row shows contours to the left and meridional held lines to the right. The 
second row shows contours to the left and the streamlines r sin6{dgv) to the right all in 
the meridional plane. 
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Fig. 10.— Regular and persistent dipole-dominated oscillations shown in the time series of 
the axisymmetric toroidal coefficients Gj (red) and Gq (blue) and the axisymmetric poloidal 
coefficients Hi (red) and (blue) describing the main dipolar and quadrupolar contribu¬ 
tions in the spherical harmonic expansion of the magnetic held in the case shown in hgure 


9 . 
















